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ABSTRACT 

This  paper  introduces  an  important  new  canonical  set  of  functions  for  solving 
Lanchester-type  equations  of  modern  warfare  for  combat  between  two  homogeneous  forces 
with  power  attrition-rate  coefficients  with  "no  offset."  Tabulations  of  these  func- 
tions, which  we  call  Lanchester-Clifford-Schlaf li  (or  LCS)  functions,  allow  one  to 
study  this  particular  variable-coefficient  model  almost  as  easily  and  thoroughly  as 
Lanchester's  classic  constant-coefficient  one.   The  availability  of  such  tables  is 
pointed  out.   We  show  that  our  choice  of  LCS  functions  allows  one  'to  obtain  important 
information  (in  particular,  force-annihilation  prediction)  without  having  to  spend  the 
time  and  effort  to  compute  force-level  trajectories.   Furthermore,  we  show  from  theoret- 
ical considerations  that  our  choice  is  the  best  for  this  purpose.   These  new  theoreti- 
cal considerations  apply  in  general  to  Lanchester-type  equations  of  modern  warfare  and 
provide  guidance  for  developing  other  canonical  Lanchester  functions  (i.e.  canonical 
functions  for  other  attrition-rate  coefficients).   Moreover,  our  new  LCS  functions 
provide  valuable  information  about  various  related  variable-coefficient  models.  Also, 
we  introduce  an  important  transformation  of  the  battle's  time  scale  that  not  only  many 
times  simplifies  the  force-level  equations  but  also  shows  that  relative  fire  effective- 
ness and  intensity  of  combat  are  the  only  two  weapon-system  parameters  determining 
the  course  of  such  variable-coefficient  Lanchester-type  combat. 


0.     Introduction 

[24] 
In  an  earlier  paper     we  developed  some  elements  of  a  mathematical  theory  for 

solving  variable-coefficient  Lanchester-type  equations  of  modern  warfare  for  combat 
between  two  homogeneous  forces  and  introduced  canonical  hyperbolic-like  Lanchester 
functions  for  constructing  their  solution.   Unfortunately,  with  only  these  previous 
results  one  is  limited  to  computing  (however  conveniently  it  may  be  done)  force-level 
trajectories  and  cannot  gain  a  real  understanding  of  qualitative  model  behavior  (e.g. 
force  annihilation)  without  the  excessive  labor  of  extensive  numerical  computations  (and 
only  then  for  specific  values  of  model  parameters) .   Since  our  earlier  work  important 
mathematical  discoveries  have  been  made  about  the  qualitative  behavior  of  this  combat 
model,  and  we  wish  to  show  how  these  new  results  allow  one  to  parametrically  analyze 
combat  modelled  by  power  attrition-rate  coefficients  (see  Section  1  below)  with  some- 
what the  same  facility  as  he  can«  study  F.  W.  LANCHESTER's     classic  constant-coeffi- 
cient model.   In  order  to  obtain  this  analysis  capability,  however,  one  must  redefine 

the  Lanchester-Clif ford-Schlafli  (or  LCS)  functions,  which  we  introduced  in  reference  24. 

[24] 
In  our  earlier  paper  (see  TAYLOR  and  BROWN    )  we  gave  various  examples  of 

hyperbolic-like  Lanchester  functions  (in  particular,  the  LCS  functions,  which  arise 

from  power  attrition-rate  coefficients  with  "no  offset") .   Subsequent  research  by 

[271 
TAYLOR  and  COMSTOCK   J  has  revealed,  however,  that  these  canonical  LCS  functions  must 

be  redefined  to  permit  force-annihilation  prediction  from  initial  conditions  without 
having  to  spend  the  time  and  effort  to  compute  force-level  trajectories.   It 
then  became  obvious  that  the  entire  topic. of  representing  the  solution  to  such  Lan- 
chester-type equations  in  terms  of  general  Lanchester  functions  (GLF)  should  be 
critically  reexamined.   Thus,  the  purpose  of  this  paper  is  to  present  new  general  con- 
siderations for  the  selection  of  canonical  Lanchester  functions  and  then  apply  this 
theory  to  the  special  case  of  power  attrition-rate  coefficients  with  "no  offset" 
(modelling,  for  example,  weapon  systems  with  the  same  maximum  effective  range)  to  obtain 
new  LCS  functions.   With  the  availability  of  tabulations  of  these  new  LCS  functions 


(see  Section  6  below) ,  one  can  study  this  model  almost  as  easily  and  thoroughly  as 
Lanchester 's  classic  constant -coef f icient  one. 

These  power  Lanchester  (i.e.  LCS)  functions  are  significant  not  only  because 
they  correspond  to  attrition-rate  coefficients  modelling  a  large  class  of  combat 
situations  of  interest  but  also  because  they  yield  valuable  information  about  other 
related  canonical  Lanchester  functions,  e.g.  the  offset  power  Lanchester  functions 
(see  Note  1  and  Section  8  below).   This  information  is,  of  course,  equivalent  to  know- 
ledge about  model  behavior  (e.g.  force  annihilation).   As  a  result  of  our  work  here 
one  can  parametrically  analyze  variable-coefficient  Lanchester-type  models  for  combat 
between  two  homogeneous  forces  with  somewhat  the  same  facility  as  the  classic  con- 
stant-coefficient one.   Such  models  are  important  for  developing  insights  into  the 
dynamics  of  combat  (see  BONDER  and  HONIg'-10',  TAYLOR^22',  and  Section  1  below). 

The  organization  of  this  paper  is  as  follows.   First,  we  present  the  variable- 
coefficient  Lanchester-type  model  that  we  study  in  this  paper.   Next,  we  discuss  the 
representation  of  the  time  history  of  the  force  levels  for  this  model  in  terms  of 
general  Lanchester  functions  (GLF) .  We  show  that  there  are  essentially  only  two  kinds 
of  GLF,  (I)  exponential-like  GLF  and  (II)  hyperbolic-like  GLF,  and  that  the  former  (I) 
provide  essential  force-annihilation-prediction  information  about  the  latter  (II). 
Then  we  explain  why  we  have  chosen  to  use  the  hyperbolic-like  GLF  to  construct  the 

model's  solution  and  why  the  power  Lanchester  (or  LCS)  functions  introduced  by  Taylor 

[24] 
and  Brown     must  be  redefined.   Next,  we  show  how  the  analysis  of,  for  example,  the 

X  force-level  equation  is  simplified  by  transforming  the  independent  variable  t   to 

normalize  the  battle's  time  scale  by  the  intensity  of  combat.   We  then  introduce  our 

new  definition  of  Lanchester-Clifford-Schlaf li  functions  and  show  how  they  arise  in 

solving  the  transformed  X  force-level  equation.   Availability  of  tabulations  of 

these  new  LCS  functions  is  discussed,  and  some  uses  of  the  tabulations  are  illustrated. 

Finally,  insights  gained  into  the  dynamics  of  combat  between  two  homogeneous  forces 

from  these  developments  are  discussed. 


1.    Variable-Coefficient  Lanchester-Type  Equations  of  Modern  Warfare 

In  this  paper  we  consider  the  following  idealized  model  for  combat  between  two 
homogeneous  forces  (see  Note  2) 

dx/dt  =  -a(t)y,  dy/dt  =  -b(t)x,     (1) 

with  initial  conditions 

x(t=0)  =  xQ,  and  y(t=0)  =  yQ, 

where   t=0  denotes  the  time  at  which  the  battle  begins,   x(t)   and  y(t)   denote  the 
numbers  of  X  and  Y  at  time  t,   and  a(t).  and  b(t)   denote  time-dependent  Lan- 
chester  attrition-rate  coefficients.   We  will  refer  to  (1)  as  variable-coefficient 
Lanchester-type  equations  of  modern  warfare  in  honor  of  the  pioneering  military  opera- 
tions research  work  of  F.  W.  Lanchester1"  '    (see  TAYLOR1  '    and  Taylor  and  Brown L   J). 
Other  forms  of  Lanchester-type  equations  appear  in  the  literature,  but  we  will  not  con- 
sider these  here  (see  DOLANSKY     and  Taylor    ).   The  Lanchester-type  equations  (1) 
yield  the  X  force-level  equation 

d2x/dt2-  {d  £na(t)/dt}dx/dt-a(t)b(t)x  -   0,  (2) 

with  initial  conditions 

x(t=0)  =  xQ,  and        { [l/a(t)]dx/dt}t=0  -  -yQ. 

Although  combat  between  two  military  forces  is  a  complex  random  process,  such 

an  idealized  deterministic  model  of  the  combat  attrition  process  is  frequently  employed 

[9] 
to  provide  insights  into  the  dynamics  of  combat  (see,  for  example,  BONDER  and  FARRELL   , 

Bonder  and  Honig^10^,  TAYLOR  and  PARRY^28\  or  WEISS^29b-   The  reader  may  consider  (1) 
to  model  combat  in  which  both  sides  use  aimed  fire  and  target  acquisition  times  are  inde- 
pendent of  the  numbers  of  f irers  and  targets  (see  Note  3) .   New  operations  research 

roi  r 121 

techniques  (see,  for  example,  Bonder  and  Farrell1   ,  and  CLARK1   J)  for  forecasting 

temporal  variations  in  fire  effectiveness  (caused  by,  for  example,  changes  in  force 

separation,  combatant  postures,  target  acquisition  rates,  etc.)  have  generated  interest 

in  such  variable-coefficient  combat  formulations. 

3 


Without  loss  of  generality,  we  may  take  a(t)   and  b(t)   to  be  of  the  form 

a(t)  =  kag(t),  and  b(t)  =  k^t),      (3) 

where  g(t)   and  h(t)   denote  the  time-varying  factors  of  a(t)   and  b(t)   such  that 
a(t)/b(t)  =  k  /k,   for  g(t)  =  h(t).   In  other  words,   k   and  k,   denote  "scale" 
factors  chosen  so  that  the  case  of  constant  coefficients  corresponds  to  g(t)  =  h(t)  = 

A  large  class  of  tactical  situations  of  interest  can  be  modelled  with  the  follow 
ing  general  power  attrition-rate  coefficients 

a(t)  =  ka(t+Ks)y,         and.  b(t)  =  kb(t+Ks+KQ)V,     (4) 

where  K  ,Kn  ^  0.   The  modelling  roles  of  K   and  K0  are  discussed  in  Taylor  and 

F241 
Brown    .   We  will  call  K   the  starting  parameter,  since  it  allows  us  to  model 

(with  u,v  £  0)  battles  which  begin  within  the  maximum  effective  ranges  of  the  two 

systems.   We  will  call  Kn  the  offset  parameter,  since  it  allows  us  to  model  (again, 

with  u,v  £  0)  battles  between  weapon  systems  with  different  effective  ranges.   Restric 

tions  that  must  be  placed  in  u   and  v,   which  are  not  necessarily  integers,  are 

discussed  below. 

Let  us  take  a  few  moments  to  motivate  our  above  notation  and  further  indicate 

[4  5] 
possible  applications  of  our  results.   Consider  BONDER'S   *   constant-speed  attack 

on  a  static  definsive  position  modelled  by 

dx/dt  =  -a(r)y  =  -a0(l-r/Ra)yy ,     dy/dt  =  -6(r)x  =  -60(l-r/R3)Vx,    (5) 

where  u,v  £  0  and  R   denotes  the  maximum  effective  range  of  the  Y  weapon  system. 

a 

Then  the  starting  parameter  and  the  offset  parameter  are  given  by 

Kg  =  (Ra-RQ)/v,  and  KQ  =  (R^-R^/v,       (6) 

where  R   denotes  the  battle's  opening  range  and  v  >  0  denotes  the  constant  attack 
speed.   Hence,   Kn,Kc  ^   0  «»  R_  £  R  ^  R  .   By  considering  (6)  and  Figure  1,  the  reader 

U     b  P         Ct         (J 

should  have  no  trouble  in  understanding  our  terminology  for  K   and  KQ.   In  the  model 
(5)  u,   for  example,  is  used  to  model  the  range  dependence  of  Y's  attrition-rate 
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coefficient  (see  Figure  2).   Observing  that  range  is  related  to  time  by  r(t)  =  R_-vt, 

we  readily  see  that  the  longest  the  battle  can  last  is  given  by  t    =  R~/v,   at  which 

max    0 

time  zero  force  separation  is  reached. 

When  the  offset  parameter  is  equal  to  zero  (i.e.   K~  =  0),  then  the  coefficients 
(4)  reduce  to 

a(t)  =  ka(t+Ks)P,         and  b(t)  =  kfe(t+Ks)V.      (7) 

We  will  refer  to  (7)  as  power  attrition-rate  coefficients  with  "no  offset."  The  purpose 

[24] 
of  this  paper  is  to  extend  our  previous  results     and  introduce  new  power  Lanchester 

functions  that  allow  more  information  to  be  more  conveniently  extracted  from  the  model 

(1)  with  coefficients  (7).   Specifically,  one  would  want  to  obtain  information  such  as: 

(Ql)   Who  will  "win"?   Be  annihilated? 

(Q2)   How  do  force  levels  decrease  over  time  and  how  many  survivors  will  the 
winner  have? 

(Q3)   How  do  changes  in  the  initial  force  levels  and  weapon  system  parameters 
affect  the  outcome?   Is  concentration  of  forces  a  good  tactic? 

(Q4)  How  long  will  the  battle  last? 

To  conveniently  answer  questions  (Ql) ,  (Q3)  ,  and  (Q4)  a  redefinition  of  the  Lanchester- 

Clifford-Schlaf li  (or  LCS)  functions  is  required.   Moreover,  not  only  are  results  for 

the  coefficients  (7)  of  interest  in  their  own  right  but  they  also  provide  much  valuable 

information  about  the  general  case  (4). 

2.     Representation  of  Solution  in  Terms  of  General  Lanchester  Functions 

In  this  section  we  discuss  how  to  construct  the  solution  to  the  X  force-level 
equation  (2)  in  terms  of  certain  basic  building  blocks  that  we  have  chosen  to  call 
general  Lanchester  functions  (GLF) .   We  feel  that  these  GLF  should  be  chosen  according 

to  the  guidelines  shown  in  Table  I.   Special  cases  of  these  general  considerations  have 

ro/i  ro"7i 

been  given  by  Taylor  and  Brown     and  Taylor  and  Comstock    . 

T  T 

Let  us  introduce  the  GLF  x  =  (x,  x„)   and  y  =  (y.  y0)  which  satisfy 

***      12        **      12 

x  =  ka(t)L£,  £  =  (l/k)b(t)Lx,      (8) 
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TABLE  I.   Requirements  for  General  Lanchester  Functions 


(Rl) :   They  can  be  used  to  construct  the  solutions  to  the  X  and  Y 
force-level  equations. 

(R2) :   They  should  be  as  "simple"  as  possible. 

(R3) :   A  given  set  of  functions  should  apply  to  as  large  a  class  of 
battles  as  possible. 

(RA) :   They  should  be  nonnegative. 

(R5) :   They  should  reduce  to  elementary  transcendental  functions  in 
special  cases  such  as  a  constant  ratio  of  attrition-rate 
coefficients. 

(R6) :   They  should  provide  as  much  information  as  possible  about  model 
behavior. 


where  without  further  specification  k  may  be  any  positive  constant  and  L  is  any 
2*2  square  matrix  such  that  L2  =  I.   [The  initial  conditions  for  (8)  are  to  be 
chosen  so  that  the  requirements  of  Table  I  are  met  and  are  discussed  below.]   Then  x 
satisfies  the  vector  equation 

x-{din  a(t)/dt}x-  a(t)b(t)x  =  0.  (9) 

In  other  words,  both  x.   and  x~   satisfy  the  X-  force-level  equation  (2),  and 
similarly  y   and  y_   satisfy  the  Y   force-level  equation. 

Let  us  now  investigate  all  the  possible  forms  for  the  above  GLF.   Intuitively, 
we  would  expect  two  possibilities  (keeping  the  requirements  of  Table  I  in  mind) : 
exponential-like  functions  (one  strictly  increasing  and  the  other  strictly  decreasing) 
and  hyperbolic-like  functions.   [We  are  reminded  of  these  two  possibilities  by  the 

well-known  constant-coefficient  results.]   Two  such  types  of  GLF  appear  in  Taylor  and 

r  9A l  r  97 1 

Brown     (hyperbolic-like  functions)  and  in  Taylor  and  Comstock     (exponential-like 

functions) .   We  will  show  that  these  are  the  only  two  possibilities  if  the  requirements 

of  Table  I  are  to  met  and  show  the  relationship  between  these  two  types  of  GLF. 

It  is  readily  shown  that  any   2x2  matrix  such  that  L2  =  I  must  take  one  of 

five  forms. 

LEMMA  1:   If  L2  =  I,   then  L  must  take  one  of  the  following  five  forms:   (A) 

(  )  with  6*0,   (B)  (       ),   (C)   (       ),   (D)   L  =  I,   or 

\[l-a2]/B   -a/  \Y   -1/        V  Y   1/ 

(E)   L  =  -I,  where  a,   3,   and  y     are  unrestricted  with  the  exception  that 
6*0. 

T 

If  in  addition  L  =  L   and   |l|  =  -1,   then  L  is  an  orthogonal  matrix  and  must  be 


of  the  form 

(cos   A  sin  <{>\ 

).  (10) 

sin  <f>        -cos  (J)/ 

T 
It  seems  reasonable  to  give  our  Lanchester  functions  symmetry  by  requiring  that  L  -  L  . 

[241 
We  observe  that  the  hyperbolic-like  GLF  of  Taylor  and  Brown     correspond  to 


L.  _  I     i.   In  order  that  (R4)  of  Table  I  be  met  for  a  constant  ratio  of  attrition- 
^"\1  0/ 


rate  coefficients,  i.e. 


a(t)  =  k  h(t), 
a 


and 


b(t)  =  k^hCt),         (11) 


[24] 


we  know  (see  Taylor  and  Brown    )  that  we  should  choose  k  =  »k,  /k   so  that 


x  =  \/ka   a(t)LjL, 


£  =  J\l\   b(t)Lx.       (12) 


The  general  solution  to  the  X  force-level  equation  (2)  is  given  by 


x(t)  =  C  x. 


(13) 


Introducing  the  2  *  2  matrix 


S(t)  = 


T 
J  (t) 


V(t)LT/  ' 


we  obtain  from  the  initial  conditions  to  (2)  that 


C  =  S 


Observing  that 


£  Iswl  - 


•T 
X  (t) 


T 
xT(t) 


XT(t)LT|    |£T(t 


(14) 


(15) 


we  readily  see  that   |s(t)|  =  constant  Vte(tn,+°°),  where   |s 
of  the  square  matrix  S.   Thus,  we  may  take 

|s(t)|  =  constant. 


denotes  the  determinant 


(16) 


Let  us  also  observe  that 


\/\  a(t)|s(t)|   =  W(X;L,x2) 


(17) 


where  W(x..  ,x_)   denotes  the  Wronskian  of  x1   and  x„. 


Now  let  us  subject  the  fundamental  system  of  solutions  x  to  linear  trans- 
formation 

x  =  Ax,  (18) 

such  that  the  form  of  the  equations  (12)  remains  invariant,  i.e. 


10 


A=    Vka  a(t)£X»  £-    ^T7^"b(t)£x,  (19) 

rtiere     L     again  is   such   that     L     =1.      If     L     is   given,    it   follows   that 


A  <* 


£  =  LAl£.  (20) 

?urthermore 

ST(t)  =  AST(t),  (21) 
;o  that 

|s(t)|  =  |A|-|s(t)|  =  constant.  (22) 

Je  also  observe  that  W(x- ,x„)  =  | A|w(x1  ,x„)  .%  Considering  the  quotient  of  the  two 


general  Lanchester  X-functions  (GLXF) 


n(t)  =  Xl/x2>  (23) 

je   see  that  under  the  linear  transformation  (18)  we  have 

dn/dt  =  {|A|/(a12n+a22)2}dn/dt.  (24) 

We  now  show  that  the  only  possible  GLF  that  satisfy  the  requirements  of  Table  I 

T 
.corresponding  to  L  =  L  )   are  the  exponential  ones  shown  in  Table  II  and  the  hyper- 
bolic ones  shown  in  Table  III.   According  to  Lemma  1  if  L2  =  I ,   then  L  must  take 
me  of  five  forms.   It  is  impossible  to  have  L  =  -I   and  satisfy  (R4)  of  Table  I  (see 
lote  4).   If  L  =  I  and  we  try  to  specify  the  "simplest"  initial  conditions  [i.e. 
ipecify  initial  conditions  such  that  the  GLF  take  the  "simplest"  form  (satisfy  (R2)  of 
'able  I)],  we  find  that  we  may  take  L  to  have  one  of  the  three  remaining  forms  (see 
fote  5).   If  we  require  that  L  be  symmetric  for  simplicity  [requirement  (R2)],  then 
.  is  an  orthogonal  matrix  with  the  form  (10).   If  in  (10)  we  take  cos  <j>  and  sin  <f 
:qual  to  -1,  0,  or  1   in  order  that  the  GLF  take  the  "simplest"  form,  we  find  that 
:he  only  two  distinct  possibilities  for  L  are  L  and  L^  as  given  in  Tables  II 
md  III  (see  argument  given  in  Note  5).   Thus,  we  have  shown  that  if  we  wish  to  con- 
struct the  solution  (13)  to  (2)  by  using  GLF  with  the  properties  given  in  Table  I, 
:here  are  essentially  only  two  possibilities:   the  exponential-like  GLF  introduced  by 

11 


Table  II.   Exponential-Like  General  Lanchester  Functions 


1.   x£(t)  =  (E+(t;Q*)   E~(t;Q*")),  ^(t)  =  (E+(t;Q*)   E~(t;Q*)) 


2«   ^<t=t0}  =  -(1/Q*     1}*    ^(t=t0}  =  (1   Q*> 


1    0' 
3.  L. 


"(        ) 

\0   -1' 


4.  |s_(t)|  =  -2 


12 


Table  III.   Hyperbolic-Like  General  Lanchester  Functions 


1.   x£(t)  =  (Sx(t)   Cx(t)),   -^(t)  =  (SY(t)   Cy(t)) 


2.     ^(t=t0)  =  (0  1),   y^(t=tQ)  =  (0  1) 


;0   1 
3.  L. 


JH  =  (i  o) 


A.  |SR(t)|  =  -1 


13 


Taylor  and  Comstock     and  the  hyperbolic-like  GLF  introduced  by  Taylor  and  Brown    . 

We  call  the  quantity  Q*,  which  appears  in  the  initial  conditions  for  the 
exponential-like  GLF,  the  parity  condition  parameter.   It  is  chosen  so  that  (see  Note  6) 

E~(t;Q*),  E^(t;Q*)  >  0  for  all  finite  t  £  tQ.  (25) 

It  may  be  considered  to  be  the  enemy  force  equivalent  of  a  friendly  X  force  of  unit 

[27] 
strength.   Taylor  and  Comstock     show  how  knowledge  of  the  parity  condition  parameter 

allows  one  to  predict  force  annihilation  from  initial  conditions  without  explicitly 

computing  force-level  trajectories.   We  observe  that  the  exponential-like  GLF  cannot 

[271 
be  computed  until  one  has  solved  the  associated  auxiliary  parity-condition  problem 

(i.e.  knows  how  to  predict  force-annihilation).   For  this  reason  and  others  (see  Taylor 

[27] 
and  Comstock    ),  the  exponential-like  GLF  are  mainly  of  theoretical  importance. 

Moreover,  in  the  next  paragraph  we  show  how  the  exponential-like  GLF  provide  valuable 

force-annihilation  information  about  the  hyperbolic-like  GLF. 

We  now  show  that  the  limiting  value  of  the  quotient  of  the  two  hyperbolic-like 

GLXF,   nH  =  S X/C„,   is  equal  to  the  reciprocal  of  the  parity  condition  parameter, 

i.e.  (30)  holds.   We  know  that  the  two  types  of  GLF  are  related  by  a  linear  transformatioi 


From  (21)  at   t  =  0,  we  have 


2«  =  ^Se*  (26) 


'0  1\     /1/Q*    1 


4  0'     W    -Q*' 


whence 


1/2    -1/(2Q*)\ 
A=[  ).  (27) 


/1/2    -1/(2Q*)\ 
Vq*/2     1/2   / 


Considering  (23)  and  (26) ,  we  see  that 

nHU)  =  {ailnE(t;Q*)+a12}/{a21nE(t;Q*)+a22}.  (28) 

Recalling  that   ri  (t;Q*)  =  E*(t;Q*)/E3(t;Q*)   and  lim  E^(t;Q*)  -  0,  we  see  that 
E         X        X  t-H«  X 
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lim  nP(t;Q*)  =  -H»,  (29) 

t->4~  fi 

so  that  by  (27)  and  (28) 

lim  n  (t)  =  1/Q*.  (30) 

t-H~  H 

rhis  result  (30)  is  highly  significant  because  it  leads  to  a  computation  method  for 
determining  Q*.   Moreover,  in  the  future  we  will  show  how  the  LCS  functions  introduced 
in  this  paper  (see  Section  5  below)  play  a  crucial  role  in  such  numerical  determinations 
Let  us  further  note  that  (24)  reads 

dnR/dt  =  {2/(Q*rig+l)2}dnE/dt,  (31) 

so  that  both  n^Ct)  and  n  (t)  are  strictly  increasing  functions  of  t,   since 
E  ti 

3nE/dt  =  2A~7ir  a(t)/{E~(t;Q*)}2. 

3,     Considerations  for  the  Choice  of  General  Lanchester  Functions 

In  Table  I  we  give  the  general  requirements  that  we  feel  should  be  placed  upon 
3LF.   These  requirements  are  motivated  by  the  properties  possessed  by  the  functions 
(namely,  the  exponential  and  hyperbolic  ones)  that  one  uses  to  construct  the  solution 
to  (1)  in  the  constant-coefficient  case. 

We  specify  (R3)  so  that  as  few  tabulations  of  GLF  as  possible  will  be  required. 

X  Y 
Consequently,  we  specify  the  initial  conditions  for  the  GLF  at   t_  =  max(tn,t_), 

X 
fhere   t_   denotes  the  largest  finite  singular  point  on  the  t-axis  for  the  X  force- 
Level  equation  (2)  (see  reference  24).   Thus,   a(t)   and  b(t)   are  positive  continuous 
functions   Vte(t0,+°°).   Since  at  most  one  of  x(t)   and  y(t)   can  vanish  in   [0,+°°) 
[see  Note  7),  we  have  what  the  mathematician  calls  a  nonoscilliatory  solution  to  (1). 
in  this  case  we  can  construct  the  solution  to  (2)  out  of  nonnegative  components  and 
rill  find  it  convenient  to  do  so. 

As  we  have  shown  above,  there  are  essentially  only  two  types  of  GLF  that 
satisfy  the  requirements  of  Table  I:   the  exponential-like  GLF  of  Table  II  and  the 
lyperbolic-like  GLF  of  Table  III.   We  feel,  however,  that  the  hyperbolic-like  functions 
ire  to  be  preferred  for  two  reasons:   (1)  they  apparently  are  more  convenient  for 
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parametric  studies  in  which  one  might,  for  example,  want  to  vary  initial  force  levels 

[24] 
or  some  measure  of  relative  fire  effectiveness    ;  and  (2)  accurate  values  of  the 

exponential-like  GLF  are,  in  general,  difficult  (in  fact,  essentially  impossible  for 

large  values  of   t)   to  determine,  since  their  initial  conditions  depend  on  the  parity 

condition  parameter  Q*   (see  Note  8).   In  terms  of  the  hyperbolic-like  GLF,  the 

[24] 
solution  to  (2)  is  given  by 

x(t)  =  x0{CY(t=0)Cx(t)-SY(t=0)Sx(t)}-y0vk^7k^"  {Cx(t-0)Sx(t)-Sx(t=0)Cx(t)  }.     (32) 

We  observe  that  for  tQ  <  0,   for  example,   C  (t=0)  >  1  and  S  (t=0)  >  0  so  that 

except  for  the  quasi-autonomous  case  in  which  a(t)/b(t)  =  constant  (see  Note  9),  the 

[24] 
solution  (32)  only  simplifies  when  tn  =  0   (see  Theorem  1  of  Taylor  and  Brown    ) 

(see  Note  10) . 

Unfortunately,  the  power  Lanchester  (or  LCS)  functions  introduced  by  Taylor  and 

[24] 

Brown     were  inappropriately  defined  to  yield  all  the  information  sought  about  the 

combat  model  (1)  with  power  attrition-rate  coefficients  (7).  In  particular,  the  time 
at  which  a  side  will  be  annihilated  cannot  be  determined  (without  the  explicit  calcu- 
lation of  the  entire  force-level  trajectories)  from  the  initial  conditions.   Subse- 

[2 
quent  work  by  Taylor  and  Comstock  has  yielded  a  theory  of  force-annihilation  prediction 

The  purpose  of  the  paper  at  hand  is  to  redefine  the  power  Lanchester  functions  in  light 
of  these  subsequent  results.   We  also  thought  it  important  to  present  the  general 
considerations  behind  this  selection  of  canonical  Lanchester  functions. 

Moreover,  the  form  of  the  LCS  functions  is  simplified  and  insight  gained  into 
the  dynamics  of  combat  by  transforming  the  battle's  time  scale.   Thus,  certain  trans- 
formations of  variables  may  be  desirable  in  the  development  of  hyperbolic-like  GLF, 
and  the  specifications  of  Table  III  should  be  interpreted  as  being  "symbolic"  and  not 
taken  literally. 
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t 
...  ds  denote  an  indefinite  integral,  denote  the  relative  effectiveness 


4.     A  Transformation  to  Normalize  the  Battle's  Time  Scale 
by  the  Intensity  of  Combat 

Let 
as  R(t) ,   i.e. 

R(t)  =  b(t)/a(t),  (33) 

and  let  K  be  an  arbitrary  constant  to  be  conveniently  chosen.   Then  Theorem  2  of 

[24] 
Taylor  and  Brown     may  be  stated  as  follows. 

THEOREM  1:   A  necessary  and  sufficient  condition  to  be  able  to  transform 
the  X  force-level  equation  (2)  by  a  transformation  of  the  independent 
variable   t   into  a  linear  second  order  ordinary  differential  equation  with 
constant  coefficients  is  that 


[dUn  R(t)}/dt]//a(t)b(t)  =  CONSTANT, 


and  then  the  desired  substitution  is  given  by 

rt 


T  =  K 


/a(s)b(s)  ds. 


We  observe  that  Theorem  1  says  that  we  can  transform  the  X  force-level  equation  to 
a  constant  coefficient  one  if  and  only  if  diHn  R(t)}/dx  =  CONSTANT.  We  also  assume 
that  the  following  condition  holds. 


CONDITION  (A) : 


r 


a(s)ds  and     b(s)ds  are  bounded  for  all  finite 
fc0 


tx  t0. 

If  Condition  (A)  is  to  hold,  then  for  the  power  attrition-rate  coefficients  (7)  we 
must  have  u,v  >  -1. 

Motivated  by  both  Theorem  1  and  the  well-known  constant-coefficient  results, 
we  introduce  the  new  independent  variable  t  defined  by 


/a(s)b(s)  ds.  (34) 


T 

t 


0 
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f  1 
By  Condition  (A)  and  the  Cauchy-Schwarz  inequality  for  integrals  (see  p.  123  of  BELLMAN1 

the  integral  in  (34)  is  well  defined  (i.e.  bounded).   The  transformation  (34)  has  an 

inverse  t(x),   since  dx/dt  >  0  V  t  >  t_.  We  also  define 

TQ  =  T(t=0).  (35) 

We  observe  that  for  tn  £  0  we  have  t_  ^  0.  Recalling  the  constant-coefficient  results 


we  will  call  the  quantity   /a(t)b(t)   the  "intensity  of  combat"  (see  also  Taylor  and 

r  28 1 

Parry    )  ;  since  the  larger  it  is,  the  more  quickly  the  battle  is  moving  towards  ter- 
mination.  The  average  intensity  of  combat  is  given  by  /a(t)b*(t)  =  (1/t)   /a(s)b(s)  ds . 


Then  we  have 


T-t0={(l/t)   /a(s)b(s)  ds}t  =  /a(t)b(t)  t.  (36) 


The  substitution  (34)  transforms  (2)  into 

d2x/dx2  +  (l/2){d  An  R(t)/dx}dx/dT  -  x  =  0,  (37) 

with  initial  conditions 

1/2 
x(t=t0)  =  xQ,     and     {R   (t)dx/dT>T=T  =  -yQ. 

•  •         u 

Theorem  1  tells  us  that  unless  (37)  is  a  constant-coefficient  equation,  it  is  impossible 

to  transform  the  X  force-level  equation  (2)  into  a  constant-coefficient  equation  by 
a  transformation  of  the  independent  variable  alone..  Also,  equation  (37)  is  highly 
significant  because  it  clearly  shows  us  that  the  course  of  combat  depends  on  just  two 
weapon-system  parameters:   (1)  R(t)  =  b(t)/a(t),   the  relative  fire  effectiveness 
(X  to  Y)  of  the  two  combatants,  and  (2)   I(t)  =  /a(t)b(t)#,   the  intensity  of  combat 
(through  equation  (34),  which  relates  I(t)   to  x) .  Both  these  parameters  may  vary 
over  time.   In  particular,  from  (37)  we  see  that  the  nature  of  temporal  variations  in 
relative  fire  effectiveness  will  have  a  significant  effect  upon  the  course  of  combat. 
For  the  power  attrition-rate  coefficients  with  no  offset  (7) ,  the  transformed 
X  force-level  equation  (37)  becomes 
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d2x/dT2+  {(2q-l)/T}dx/dT-x  =  0,  (38) 

with  initial  conditions 

x(t=tq)  =  xQ,   and   { (T/2)2q_1dx/dT}T=T  =  -yQvir^(>/k^/(u+v+2)  )2q"\ 

where  q  =  (v+l)/()j+v+2)   and 

t  =  x(t)  =  (2/kl^/(u+v+2))(t+Ks)(y+V+2)/2.  (39) 

Hence,   xQ  =  (2/k  lc/(u+v+2))Kg   V     .   Let  us  observe  that  Vu,v  >  -1  we  have 
0  <  q  <  1.   Furthermore,   q  >  1/2  *  dR/dt  >  0,   i.e.   R(t)   is  a  strictly  increasing 
function  of  time. 

5.     Lanchester-Clifford-Schlafli  Functions 

Consider  the  function  F  (x)   defined  by  the  power  series 

00 

F  (x)  =  r(o)  I      (x/2)2k/(k!r(k-Hx)}.  (40) 

k=0 

For  a  ^  0,-1,-2,...   the  radius  of  convergence  for  F  (x)   is  infinite  by  the  ratio 
test  for  convergence  of  power  series  (see,  for  example,  KN0PP    ).   Hence,   F  (z) 
is  an  entire  function  of  the  complex  variable   z  =  x+iy  with  an  essential  singularity 
at  the  point  at  infinity.   Now  consider  the  function  H  (x)   defined  by  the  infinite 


series 


Observing  that 


H  (x)  =  r(o)  I      (x/2)2(k+a)/{k!r(k-kx+l)}.  (41) 

k=0 


Ha(x)  =  (l/a)(x/2)2aFa+1(x),  (42) 


we  see  that  for  a  >  0  the  infinite  series  (41)  is  uniformly  convergent  on  compact 
subsets  of  the  complex  plane.   From  (42)  we  can  readily  deduce  the  recursive  relation 

F0(x)  =  Fa+1(x)  +  {(x/2)2/[a(a+l)]}Fa+2(x).  (A3) 

We  will  call  the  functions  F  (x)   and  H  (x)   Lanchester-Clifford-Schlafli  functions 

a  a     

(see  Note  11).   Other  properties  are  readily  deduced  and  are  given  in  Table  IV. 
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Table  IV.   Properties  of  the  LCS  Functions  F  (x)   and  H  (x) 


1-?a 
1.     dF  /dx  =  (x/2)    H  (x) 


2.     dHa/dx  =  (x/2)2a  ^(x) 


3.     F  (x)F1   (x)-H  (x)tf.   (x)  =  1     Vx 

av  '    1-a  '        av   l-av  ' 

where  a  is  not  an  integer  (including  zero) 


4.     F  (x=0)  =  1 


5.     H  (x=0)  =0     for   a  >  0 


6.     dF  /dx(x=0)  =  0 

a 


7.     {(x/2)1'2adH  /dx}  .  =  1 

a    x=0 


8.     Fi/2^  =  cosn  x 


9.     H  y7(x)  =  sinh  x 
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The  function  F  (x)   satisfies  the  second  order  ordinary  differential  equation 

d2F  /dx2+{(2a-l)/x}dF  /dx-F  =0,  (44) 

with  initial  conditions 

Fa(x=0)  =1,  and  dF  /dx(x=0)  =  0, 


while  H  (x)   satisfies 
a 


d2H  /dx2  -  {(2a-l)/x}dH  /dx  -  H  =  0,  (45) 

a  a      a 

with  initial  conditions  (for  a  >  0) 

Ha(x=0)  =  0,  and       { (x/2)1-2otdH  /dx}  =Q  =  1. 

Thus,   {F  ,H-   }   is  a  fundamental  system  of  solutions  to 

a     1-a  J 

d2F/dx2  + {(2a-l)/x}dF/dx-F  =  0,  (46) 

1  —  On 

with  Wronskian     W(F    ,H,      )   =    (x/2)  .      Let  us   observe  that    (see  Table  III) 

Cy(t)    =  F(T(t)),        S„(t)   =   {/k~kT/(M+v+2)}2q~1H    (x(t)),  (47) 

X  q  X  at)  p 

Cy(t)   =  Fp(x(t)),        Sy(t)   =   {/irk^/(y+v+2)}1"2qHq(T(t)),  (48) 
where     p  =   1-q.      If  we  define 

To(x)   =  H1_a(x)/Fa(x),  (49) 
then 

nR(t)  =  Tx(t)  =  Sx(t)/Cx(t)  =  {v^k^/(y+v+2)}2q"1Hp(T(t))/Fq(T(t)),  (50) 

where  T  (t)   denotes  a  hyperbolic-like  GLF.   Observing  that   lim  x(t)  =  +°°,  we  see 

A. 

t-*+°° 

that  T  (x)   is  a  strictly  increasing  function  of  x  Vxe[0,4-»)   and 


a 


with 


0  £.   T  (x)  <  r(l-a)/r(a)     for   0  *  x  <  -+<»,  (51) 

a 


lim  TQ(x)  =  r(l-a)/r(a),  (52) 

x-*-Ho 
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[271 
since  by  the  value  of  Q*  determined  by  Taylor  and  Comstock     for  the  power  attrition- 
rate  coefficients  (7),  denoted  as  Q*(y,v,K  =0) ,  we  have  (see  (30)  and  (50)) 

lim  Tx(t)  =  l/Q*(vi,v,Ko=0)  =  {/k^/(u+v+2)}2q~1r(p)/r(q). 
t-H-» 

Comparing  (38)  and  (46),  we  see  that  the  solution  to  (38)  is  given  by 

x(t)  =  x_{F  (t  )F  (x(t))-H  (t  )H  (x(t))} 
U   p   0   q         q   0   p 

"y0^a7k^{/k^/(u+\H^^  (53) 

The  time  to  annihilate  X  is  determined  by  x(t=tY)  =  0  and  thus 

A 


T  (x(4»  -  {x0F  (T0)+y0/T7k;(AX/(^v+2))<'-PH  (t0)}/ 


qx  *  X"      0  p^  0'  'U  a   b  '   a  b  "  "  p 

{xOHq(T0)+>'0,^71S(,^A/(l,+V+2))q"PFq(T0)}'      (54) 
where  T  (t)   is  given  by  (49)  and  from  (51) 

0  *  T  (t)  <  T(p)/r(q).  (55) 

For  tQ  =  -Kg  =  0,   (54)  simplifies  to 

Tq(x(t*))  =  (x0/y0)y^7ka"(v^k^'/(y+v+2))p"q.  (56) 

From  (54)  and  (55)  we  may  deduce  the  following  theorem: 

THEOREM  2:   Consider  combat  between  two  homogeneous  forces  described  by 
(1)  with  power  attrition-rate  coefficients  (7) .   Assume  that  these 
equations  hold  for  all  time  and  that  Y   "wins"  when  x(tf)  =  0  with 
y(t  )  >  0.   Then  Y  will  win  if  and  only  if 

r(p){x0Hq(T0)+y0^r7k^(/k^/(u+V+2))<l-t>Fq(T0)}. 
For  t  =0   (i.e.  Kc  =  0  and  x„  =  0) ,  Y  will  win  if  and  only  if 

0  b  U 

xQr(q)  <  y0/ka/kb(/kakb/(u+v+2))q"pr(P) . 


22 


6.     Tabulations  of  LCS  Functions 

Tabulations  of  the  Lanchester-Clif ford-Schlafli  functions  are  available  in  two 
of  the  authors'  reports,  also  available  from  the  National  Technical  Information  Service 
(see  references  25  and  26) .   These  reports  contain  five-decimal-place  tables  of  the 
hyperbolic-like  LCs  functions  F  (x) ,   H    (x) ,   and  T  (x)   for  values  of  the  argu- 
ment x  =  0.00(0.01)2.00(0.1)10.0  and  various  values  of  the  order  a.   The  short 
table^25^  contains  tabulations  for  a  =  1/2,1/3,2/3,1/4,3/4,1/5,2/5,3/5,4/5,3/7,  and 
4/7  corresponding  to   p,v  =  0,1,2,3;   while  the  longer  table     contains  tabulations 
for  a  =  1/2,1/3,2/3,1/4,3/4,1/5,2/5,3/5,4/5,2/7,3/7,4/7,5/7,4/9,5/9,3/11,5/11,6/11, 
8/11,5/13,8/13,5/17,12/17,5/21,  and  16/21  corresponding  to  y,v  =  0,1/4,1/2,1,1  1/2, 
2,3.   As  we  have  seen  above  in  Section  1  (see  (5)  and  Figure  2),  such  values  of  u  and 
v  allow  one  to  analyze,  for  example,  a  wide  variety  of  range  capabilities  for  weapon 

systems  in  Bonder's   '    constant-speed  attack  model  (5).   These  tables  have  been 

[241 
calculated  by  the  recursive  methods  given  in  Section  8  of  Taylor  and  Brown 

A  representative  tabulation  of  the  hyperbolic-like  LCS  functions  F  (x) , 

H-   (x) ,   and   T  (x)   for  a  =  3/5,   similar  to  those  that  appear  in  references  25  and 

26,  is  given  in  Tables  V  and  VI.   The  values  of  the  argument  x  are  the  same  as  those 

used  for  the  tabulation  of  the  hyperbolic  functions  by  ABRAMOWITZ  and  STEGUN1   .   We 

observe  from  Table  VI  and  (52)  that  the  limiting  value  of  T  (x)   as  x  ■*■   +00   (here 

a  =  3/5)  is  quickly  reached,  with  three-decimal-place  agr.ee/ient  by  x  =  4.5. 

7.     Numerical  Examples 

In  this  section  we  examine  a  couple  of  numerical  examples  to  show  some  of  the 
insights  that  may  be  gained  into  the  dynamics  of  combat  between  two  homogeneous  forces 
from  our  new  results.   As  in  references  21  and  24,  we  consider  S.  Bonder's      model 
(5)  for  the  constant-speed  attack  against  a  static  defensive  position.   We  will  focus 
on  the  new  results  of  this  paper  [in  particular,  the  prediction  of  battle  outcome  from 
initial  conditions  without  explicitly  computing  the  force-level  trajectories  (cf. 
questions  (Ql)  and  (Q4)  of  Section  1) ] .   From  the  input  data  given  in  Table  VII,  we 
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Table  VII.   Input  Data  for  Numerical  Examples 


u  =  1,   v  =  2 


an  =  0.06  X  casualties/minute/Y  unit 


3_  =  0.6  Y  casualties/minute/X  unit 


R  =  R_  =  2000  meters 
a    3 


v  =  5  miles/hour 
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compute  the  parameter  values  shown  in  Table  VIII.  We  observe  from  Tables  VI  and  VIII 
the  predicted  agreement  between  T(l-a)/r(a)   and  the  limiting  value  of  T  (x)   as 
x  ■*■  +°°   (see  (52))  for  a  =  q  =  3/5.   We  now  consider  two  cases:   (I)   R_  =  2000 

meters,  and  (II)   Rn  =  1250  meters. 

T211 
When  RQ  =  2000  meters  (see  Figure  3  of  Taylor1    )  ,  we  have  K  =  0  and 

T.  =  0.   The  maximum  time  that  the  battle  can  last  is   t    =14.91  minutes,  since  at 
0  max 

this  time  the  attackers  reach  their  final  objective  (i.e.  the  defensive  position). 

We  now  consider  the  qualitative  behavior  of  the  u=  1,   v  =  2   force-level  trajectory 

[21] 
shown  in  Figure  3  of  Taylor    .   Theorem  2  tells  us  that  X  can  be  annihilated 

*>  x n/yn  <  0.420.   By  (56)  the  annihilation  time  of  the  X  force  is  given  by 

T  (T(tJ))  =  3.544  xrt/y_.   For  x.  =  10,   y_  =  30,  we  have  T  (t*)  =  1.18122  so  that 
q     A  UU  U  U  *1   X 

from  Table  V  (using  linear  interpolation)  we  obtain  T  =  1.009.   Hence,  (39)  yields 

A. 

t  =  14.24  minutes  and  r  =  89.8  meters.   Further  results  are  given  in  Table  IX. 

A  A 

[241 
When  R_  =  1250  meters  (see  Figure  3  of  Taylor  and  Brown    ) ,  we  have 

K0  =  5.5923  minutes,   T_  =  0.0975,   and   t    =  9.32  minutes.  Jn  this  case  (again, 
S  0  max 

for  u  =  1,   v  =  2),   X  can  be  annihilated  «  xn^yn  <  °*3^2  witn  [from  (54)]  the 
annihilation  time  of  the  X  force  given  by  T  (t*)  =  (3.565  u  +0.223)/(0.156  u  +1.004), 

Q   A  U  U 

where  u_  ■  x^^n"   Some  further  numerical  results  are  given  in  Table  X.   Again,  these 
parametric  results  should  be  contrasted  with  the  single  u  =  1,   v  =  2  force-level 
trajectory  shown  in  Figure  3  of  reference  24. 

8.    Discussion 

In  Section  7  above  we  have  seen  how  our  new  definition  of  power  Lanchester 
functions  (guided  by  the  general  requirements  for  GLF  given  in  Table  I)  allows  one  to 
conveniently  obtain  much  valuable  information  about  the  model  (1)  with  attrition-rate 
coefficients  (7)  without  explicitly  computing  the  entire  force-level  trajectories  (see 
Note  12).   Previously  we  were  limited  to  only  computing  force-level  trajectories. 
Now  we  can  tell  who  is  going  to  be  annihilated  and  when  without  explicitly  computing 
the  trajectories  (see  Note  13).   Not  only  did  we  answer  questions  about  qualitative 
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Table  VIII.   Parameter  Values  for  Numerical  Examples 


k  =  4.0233*10   X  casualties/ (minute)y/Y  unit 
a 


k,  =  2.6979x10"  Y  casualties/  (minute)  V/X  unit 


p  =  2/5,   q  =  3/5 


r(p)/r(q)  =  1.48951 


KQ=0 
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Table  IX.   Annihilation  of  the  X  Force  as  a  Function 
of  the  Initial  Force  Ratio  for  RQ  =  2000  meters 


V^O1 

t    (minutes) 

A 

r    (meters) 

0.333 

14.24 

89.8 

0.250 

11.61 

443.2 

0.200 

10.19 

633.2 
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Table  X.   Annihilation  of  the  X  Force  as  a  Function 
of  the  Initial  Force  Ratio  for  RQ  =  1250  meters 


(xQ/yQ)  tx (minutes)  rx(meters) 

0.333  10.63             t 

0.250  7.56                                 235.9 

0.200  6.17                                  422.8 


t 
t  =9.32  minutes     and     x,  =  x(r=0)   =   1.35 

max  f 


30 


model  behavior  (e.g.  force  annihilation)  for  specific  values  of,  for  example,  initial 
force  levels  but  also  for  a  range  of  values  of  the  initial  force  ratio  (i.e.  parametric 
analysis  of  model  behavior) .   The  results  of  this  paper  may  be  used  for  other  para- 
metric analyses  (see  Note  14),  e.g.  parametric  dependence  of  battle  outcome  on  attri- 

[241 
tion-rate  coefficients.   Thus,  our  extension  of  past  results     allows  one  to  develop 

important  insights  into  the  dynamics  of  combat  between  two  homogeneous  forces  with 

T  2  S  261 
temporal  variations  in  fire  effectivenesses.   With  the  availability    '    of  tabula- 
tions of  the  LCS  functions,  one  can  now  analyze  such  combat  modelled  by  the  power 
attrition-rate  coefficients  (7)  with  somewhat'  the  same  facility  as  he  can  for  the 
constant-coefficient  case  (see  Note  15)  and  thus  aid  in  parametric  analyses. 

In  his  classic  1914  paper     Lanchester  assumed  that  the  combatants'  fire 
effectivenesses  (as  expressed  by  Lanchester  attrition-rate  coefficients)  were  constant 
over  time  and  deduced  his  famous  square  law 

B{X2-x2(t)}  =  a{y2-y2(t)},  (57) 

where  a  and  8  denote  constant  attrition-rate  coefficients.   It  follows  from  (57) 
that  (provided  there  is  no  "time  limit"  for  the  battle) 


X  will  be  annihilated  »  xr/yn  <  *°^&m  ^58^ 

Thus,  we  see  that  equality  of  Lanchester-type  fighting  strengths  depends  on  two  param- 
eters:  (I)   initial  force  ratio,  and  (II)   relative  effectiveness.   When  the  timing 
of  military  actions  is  considered,  we  add  a  third  parameter,  the  intensity  of  combat  = 
ya&y      to  this  list  of  significant  combat  parameters.   No  such  simple  relationship  like 
the  square  law  (57),  which  yields  (58),  holds  in  general  for* variable  attrition-rate 
coefficients  when  a(t)/b(t)  J*  constant.   However,  by  transforming  the  independent 
variable  t   to  normalize  the  battle's  time  scale  by  the  intensity  of  combat,  we  found 
(see  equation  (37))  that  the  course  of  such  variable-coefficient  combat  depends  on 
only  two  weapon-system  parameters:   (I)   relative  fire  effectiveness,   R(t)  =  b(t)/a(t), 


and   (II)   intensity  of  combat,   I(t)  =  /a(t)b(t) .   Moreover,  we  extended  (58)  to  combat 
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modelled  with  the  power  attrition-rate  coefficients  with  "no  offset"  (7)  (see  Theorem  2) 
This  is  the  first  time  that  such  a  generalization  of  the  square  law  has  been  obtained 
for  the  variable-coefficient  Lanchester-type  model  (1)  with  a(t)/b(t)  f   constant.  We 
observe  that  for  K  >  0  this  "exact"  outcome-prediction  relation  (i.e.  necessary 
and  sufficient  condition  for  force  annihilation)  involves  higher  transcendental  functions 

(here,  the  LCS  functions)  and  is  complementary  to  the  sufficient  condition  (involving 

r  or  1 

only  elementary  functions)  given  by  Taylor  and  Parry     for  K„  >  0. 

\5   71       T121  \2   91 

Work  by  BONDER1  '   ,  Clark1   J,  and  others1  '  J  on  the  prediction  of  Lanchester 

[241 
attrition-rate  coefficients  (see  Taylor  and  Brown     for  further  discussion  and 

references)  has  generated  interest  in  variable-coefficient  Lanchester-type  models. 
Interest  in  the  power  attrition-rate  coefficients  with  "no  offset"  (7)  is  provided  by 
S.  Bonder's   '   model  (5)  and  his  examination  of  predicted  attrition-rate  for  various 
weapon  systems  (see  pp.  196-200  of  reference  9).   However  useful  our  results  may  be 
in  their  own  right,  they  have  far  greater  import:   (I)   they  are  a  model  for  the  treat- 
ment of  other  Lanchester  functions  and  their  tabulations,  and  (II)  they  may  be  used  in 

[271 
the  numerical  determination  of  the  parity-condition  parameter      Q*  for  related 

attrition-rate  coefficients  (for  example,  (4)  with  KQ  >  0) .   In  the  future  we  will 
show  how  our  tabulations  of  the  LCS  functions  play  a  key  role  in  the  numerical  deter- 
mination of  the  parity-condition  parameter  Q*  for  the  general  power  attrition-rate 

coefficients  (4)  with  positive  "offset"  (i.e.  KQ  >  0) . 

[241 
We  have  extended  our  mathematical  theory     of  variable-coefficient  Lanchester- 
type  equations  of  "modern  warfare"  for  combat  between  two  homogeneous  forces  in  order 
to  be  able  to  more  thoroughly  analyze  such  models.   The  classic  ordinary  differential 
equation  theories  (see,  for  example,  HILLE    )  were  inadequate  to  supply  all  the 
answers  sought  about  such  combat  models  (cf.  questions  (Q1)-(Q4)  in  Section  1  above). 
The  mathematical  theory  of  the  model  (1)  with  coefficients  (7)  is  now  nearly  as  complete 
as  that  of  the  constant-coefficient  model.   Such  results  as  we  have  given  here  are  very 
useful  for  understanding  the  dynamics  of  combat,  i.e.  how  the  trading  of  casualties 
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will  be  projected  over  time.   H.  K.  WEISS     has  emphasized  that  such  a  simplified 
model  of  a  combat  situation  is  particularly  valuable  when  it  leads  to  a  clearer 
understanding  of  significant  relationships  that  would  tend  to  be  obscured  in  a  more 
complex  model.   As  is  always  the  case,  however,  the  insights  geined  into  combat 
dynamics  are  no  more  valid  than  the  models  themselves. 

9.    Summary 

In  this  paper  we  have  introduced  new  mathematical  functions  (Lanchester-Clifford- 
Schlafli,  or  LCS,  functions)  that  allow  important  information  (in  particular,  force- 
annihilation  prediction)  to  be  obtained  without  explicitly  computing  force-level 
trajectories  for  the  variable-coefficient  Lanchester-type  model  (1)  with  power  attri- 
tion-rate coefficients  with  "no  offset"  (7).   Our  development  was  based  on  new  theoretical 
considerations:   we  gave  a  new  general  discussion  of  representing  the  solutions  to  the 
X  and  Y   force-level  equations  in  terms  of  general  Lanchester  functions  (GLF)  and 
gave  the  general  properties  that  these  GLF  should  possess;  we  showed  that  there  are 
essentially  only  two  kinds  of  GLF  that  satisfy  these  requirements  (exponential-like 
GLF  and  hyperbolic-like  GLF)  and  that  the  hyperbolic-like  functions  are  to  be  preferred. 
Moreover,  the  exponential-like  GLF  are  an  essential  theoretical  construct,  since  they 
play  a  key  role  in  determining  force-annihilation-prediction  conditions  (i.e.  showing 
that  the  reciprocal  of  the  parity  condition  parameter  is  equal  to  the  limiting  value 
of  the  quotient  of  two  hyperbolic-like  general  Lanchester  X-functions) .  We  stressed 
that  such  building  blocks  should  be  chosen  to  yield  as  much  information  as  possible 
about  the  model  (and  as  conveniently  as  possible).   We  saw  that  the  analysis  of,  for 
example,  the  X  force-level  equation  was  facilitated  by  transforming  the  battle's 
time  scale  and  that  the  only  two  weapon-system  parameters  affecting  the  course  of 
combat  are  the  relative  fire  effectiveness  and  the  intensity  of  combat.   These  results 
extended  and  unified  our  mathematical  theory  of  variable-coefficient  Lanchester-type 
equations  of  "modern  warfare"  (see  reference  24) . 
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We  then  applied  our  general  mathematical  theory  to  the  special  case  of  combat 
modelled  by  power  attrition-rate  coefficients  with  "no  offset."  Our  new  definition 
of  Lanchester-Clifford-Schlafli  (LCS)  functions  was  required  for  these  power  attrition- 
rate  coefficients  in  order  to  answer  questions  about  battle  outcome  without  explicitly 
computing  force-level  trajectories  (i.e.  to  predict  battle  outcome/force  annihilation). 
The  mathematical  theory  of  this  variable-coefficient  Lanchester-type  models  of  "modern 
warfare"  (modelling,  for  example,  weapon  systems  with  the  same  effective  range)  is 
now  nearly  as  complete  as  that  of  the  constant-coefficient  model.   With  tabulations 
of  the  new  LCS  functions  now  available,  one  can  study  this  variable-coefficient  model 
almost  as  easily  and  thoroughly  as  Lanchester' s  classic  constant-coefficient  model. 


Notes 

1.  Following  terminology  introduced  in  reference  24,  we  will  refer  to  Lanchester 
functions  corresponding  to  the  power  attrition-rate  coefficients  (4)  with  Kq  >  0  as 
offset  power  Lanchester  functions  (see  Section  1).   The  power  Lanchester  (i.e.  LCS) 
functions  correspond  to  K  =0,   i.e.  to  the  power  attrition-rate  coefficients  with 
"no  offset"  (7). 

2.  The  equations  (1)  are  only  valid  for  x,y  >  0.   The  first,  for  example,  becomes 
dx/dt  =  0  for  x  =  0. 

3.  Further  information  on  sets  of  circumstances  that  have  been  hypothesized  to  yield 
the  combat  equations  (1)  (with  constant  coefficients)  may  be  found  in  BRACKNEY 

and  Weiss 

4.  It  is  impossible  for  all  a(t),b(t)  >  0  satisfying  Condition  (A)  to  have 

x.,y.  >  0  for  i  -  1,2  and  all  t  >  trt  such  that  x.  =  -/k,/k  a(t)y.   and  y.  = 
11  U  i      d  a      i        i 

-A.  /k,  b(t)x.   and  x..  ,x~  are  linearly  independent.   If  it  were  possible,  then 
|s(t)j  =  -x1y2  +  x_y1  =  C-  and  without  loss  of  generality  we  may  take  C_  >  0. 
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Introducing  u.  =  x  /y   and  A  =  u  -  u- ,  we  would  have   |s(t)|  =  y-.y9A  =  C  >  0  and 
I   =  J\l\   b(t)(U;L+u2)A  >  /k^/^  b(t)A2  >  0  for  t  >  tQ  so  that  A(t)  >  0  is  strictly 
increasing  for   t  £  tQ.   It  would  follow  that  1/A(t0)  -  1/A(t)  £  /k_  /kr 


which  is  impossible  for  b(t)  such  that   lim 

t-*+°° 


t 

b(s)ds, 

■t  v0 

b(s)ds  =  +°°. 


5.  To  keep  x_,  x.,  y..  ,  and  y   as  "simple"  as  possible  we  specify  that  their  initial 
values  at   t_   to  be  either  0  or   1.   In  order  that   |s(t)|  =  x1y/?-x_y1  ^  0,  we 
must  therefore  have  either  x1(t=tQ)  =  y_(t=t0)  =  0  and  x2(t=tQ)  =  y..  (t=t0)  =  1,   or 
x1  (t-t0)  =  y„(t=t_)  =  1  and  x_(t=tn)  =  y..  (t=t-)  =  0.   We  consider  the  first  possi- 
bility with  similar  arguments  holding  for  the  second.   In  this  first  case  examination 
of  the  differential  equations  with  initial  conditions  shows  us  that  x.  =  S  ,   x  =  C  , 

J.     A     Z     X 

y  =  C  ,   and  y_  =  S  ,   i.e.  the  functions  coincide  with  the  hyperbolic-like  GLF  of 
Table  III.   Thus,  we  need  not  consider  L  =  I,   since  the  same  results  may  be  obtained 
by  using  another  one  of  Lemma  l's  feasible  forms  for  L. 

6.  We  conjecture  that  some  condition  like  lim  x(t)  =  +°°  is  sufficient  to  guarantee 

t-H-°° 
that  Q*  is  unique. 

7.  This  intuitively  obvious  result  may  be  proved  by  observing  the  identity 

{b(a)x2(o)  +  a(o)y2(a)}da  =  x(s)y(s)  -  x(t)y(t) . 
s 
A  less  obvious  fact  is  that  unless  at  least  one  of 


a(t)dt   and 
0 


b(t)dt  is 
0 


unbounded,  then  neither  x(t)   nor  y(t)   need  ever  be  annihilated  (see  Hille    ). 

As  an  example  of  this  situation,  consider  the  battle  with  attrition-rate  coefficients  (7), 

K  >  0,   and  u  ■  V  <  -1.   Then  x(t)  =  xQ  cosh  0 (t)  -  yQ/k  /1c  sinh  9  (t)  ,  where  6(t)  = 

{-l/(v+l)Hl/K  "(v+1)-l/(t+Kc)~(v+1)  .   Let   6  =  {-l/(v+l)}{l/K _"(v+1)}  >  0  and 

finite.   It  follows  that  v  <  -1  and  K  >  0  can  be  chosen  so  that  even  though 

x«  <  yrt/k  /k,   we  have  lim  x(t)  =  x_  cosh  6  -  yrt/k  /k_  sinh  9   >  0. 
0^0  a  D  ,  0      »   0  a  d      °° 
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8.  In  general,  the  value  of  Q*  will  not  be  known  exactly.   Unfortunately,  errors 

in  the  initial  conditions  for  the  exponential-like  GLF  become  exponentially  magnified 

+  +        - 

over  time.  The  situation  is  even  worse  for  n_(t;Q*)  =  E_  (t;Q*)  =  E  (t;Q*)/E  (t;Q*) , 

which  is  used  to  determine  the  time  that  the  X  force  will  be  annihilated. 

[22]  [23] 

9.  The  term  quasi-autonomous  was  coined  by  Taylor     (see  also  TAYLOR    )  to  denote 

a  system  of  differential  equations  that  may  be  transformed  to  an  autonomous  system  (see 

risi 

for  example,  p.  163  of  PETROVSKI    )  by  a  change  of  the  time  scale.   Special  cases 

[9] 
of  such  Lanchester-type  equations  have  been  considered  by,  for  example,  Farrell    and 

TAYLOR    .   More  general  (possibly  nonlinear)  quasi-autonomous  Lanchester-type  equa- 

rooi  fo^i  r  o  /  1 

tions  have  been  studied  by  Taylor         (see  also  Note  4  of  Taylor  and  Brown    ). 

10.  If  we  were  to  specify  the  initial  conditions  of  the  GLF  at  t  =  0  instead  of 

t  ■  t_,   then  (32)  would  reduce  to  x(t)  =  x_Cv(t)  -  y_/k  /k~~  sL(t)  .   However,  when  the 
U  (J  X     '  U  a  b   a 

initial  conditions  for  the  hyperbolic-like  GLF  are  not  given  at  tn,  a  separate  tabu- 

•*> 
lation  of,  for  example,  C Y(t)  must  be  used  for  each  different  value  of  tn 

(i.e.   Cx=  Cx(t;t0)). 

11.  Although  the  solution  to  the  X  force-level  equation  (2)  with  the  power  attrition- 
rate  coefficients  (7)  may  be  expressed  in  terms  of  known  higher  transcendental  functions 

roil  ro/i  ro7i 

(see  Taylor    ,  Taylor  and  Brown    ,  and  Taylor  and  Corns tock    ) ,  we  have  chosen 

to  introduce  the  LCS  functions,  since  tabulations  of  these  other  functions  are  not 

available  for  the  full  range  of  parameter  values  of  interest  in  Lanchester  combat 

theory.   For  example,  we  can  construct  such  solutions  with  modified  Bessel  functions 

of  the  first  kind  of  fractional  order,  but  tabulations  of  these  (see,  for  example, 

Abramowitz  and  Stegun   )  only  exist  for  a  restrictive  set  of  values  of  the  order  p 

(i.e.   p  -  ±1/4,  ±1/3,  ±1/2,  ±2/3,  ±3/4),  where  p  =  (u+l)/(u+v+2) .  Furthermore, 

tabulations  of  functions  corresponding  to  the  quotient  of,  for  example,  two  GLXF  do 

not  apparently  exist  at  all.   Consequently,  we  have  introduced  our  new  LCS  functions, 

which  provide  much  of  the  information  desired  about  such  battles.   The  naming  of  our 

LCS  functions  follows  from  the  facts  that  a  function  similar  to  F  (x)  was  introduced 

a 
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by  LUDWIG  SCHLAFLI1   J  (1814-1895)  in  1867,  while  a  related  one  appears  in  a  posthumous 
fragment  of  the  great  English  geometer  William  Kingdon  Clifford  (1845-1879)  (see  pp.  343- 
348  of  CLIFFORD ^13]). 

r  14 1 

12.  In  his  well-known  survey  paper  on  the  Lanchester  theory  of  combat,  Dolansky 

suggested  the  development  of  outcome  predicting  relations  without  solving  in  detail 
and/or  computing  force-level  trajectories  as  one  of  several  problems  for  future  research. 
Our  Theorem  2  is  a  step  towards  this  problem's  resolution  (see  also  references  22,  27, 
and  28). 

13.  Bonder  and  Honig     point  out,  however,  that  force  annihilation  may  not  be  the 

appropriate  criterion  for  evaluating  many  military  operations,  especially  when  force 

[91 
annihilation  does  not  occur.   See  pp.  192-242  of  Bonder  and  Farrell    for  a  detailed 

Lanchester-type  analysis  of  an  attack  situation  for  which  other  "end  of  battle  condi- 
tions" play  the  major  role  in  the  evaluation  process.   Nevertheless,  it  is  of  interest 
to  know  when  and  why  force  annihilation  will  occur. 

14.  S.  BONDER    has  suggested  that  an  increased  emphasis  be  placed  on  parametric 
analyses  in  systems  analysis  studies  (see  pp.  21-22  of  reference  8). 

15.  One  significant  exception  is  that  the  outcome  of  fixed-force-level-breakpoint 
battles  (for  example,   Y   "wins"  when  xf  =  x(tf)  =  x  p  but  yf  >  yRp,  where  t  , 
xr>  y*  denote  final  values  and  xDn  denotes  X's  breakpoint)  with  x  D  y  „  >  0  and 

t    1  Br  Br ,  Br 

a(t)/b(t)  f   constant  cannot  apparently  be  analyzed  in  the  manner  described  in  this 

[271 
paper  (see  Taylor  and  Comstock    ). 
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